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Abstract 

We consider scattering and capture of circular cosmic strings by a 
Schwarzschild black hole. Although being a priori a very simple axially sym- 
metric two-body problem, it shows all the features of chaotic scattering. In 
particular, it contains a fractal set of unstable periodic solutions; a so-called 
strange repellor. We study the different types of trajectories and obtain the 
fractal dimension of the basin-boundary separating the space of initial con- 
ditions according to the different asymptotic outcomes. We also consider the 
fractal dimension as a function of energy, and discuss the transition from order 
to chaos. 
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Chaos in general relativity and cosmology is by now a well-established subject. Some 
studies, pioneered by the work of Hawking and Page 0], have been concerned with the 
behavior of solutions to the Einstein equations themselves; the most famous example now 
being the mixmaster universe (see for instance [Q, and references given therein). In other 
studies, the subject of interest has been the behavior of test-particle trajectories in black 
hole spacetimes (see for instance In the latter case, since point-particle motion is com- 

pletely integrable in the generic Kerr- Newman background 0, it was necessary to consider 
quite complicated multi black hole spacetimes (typically of the Majumdar-Papapetrou type 
to obtain chaotic point-particle dynamics. Taking into account that these systems 
represent (at least) three-body problems, and comparing with Newtonian dynamics, it was 
certainly no surprise that chaotic dynamics appeared. (It should be mentioned, however, 
that in the two black hole case, chaos is really a relativistic effect [^.) The main interest in 
these systems therefore was also to use and further develop coordinate-invariant descriptions 
and measures of chaotic behavior @JT^ , suitable for general relativity where space and time 



are not absolute. 

In the present paper, we study a conceptually simpler and more symmetric system in 
general relativity that nevertheless, and maybe somewhat surprisingly, leads to chaotic be- 
havior. We consider the two-body problem of a circular test-string in the background of 
a Schwarzschild black hole. The circular string is taken to be coaxial with the black hole 
and is allowed to oscillate in its plane and to propagate in the direction perpendicular to its 
plane, as illustrated in Fig. |l|. The combined system is therefore axially symmetric. 

The physical picture we have in mind is that of a cosmic string nearby an astrophysical 
black hole. For a GUT string with string tension of the order G/i ~ 10~^ |]1T|, the test-string 
approximation should be valid even for a string initially up to (say) 10^ times longer than 
the black hole horizon radius. Moreover, for a GUT string we can use the leading order 
thickness approximation, where the cosmic string is described by the Nambu-Goto action 



S = —jj. J dr da \J — det [Gab\ , (1) 

where Gah = gfiuX^a-^t is the induced metric on the string worldsheet. In this case, the 
string equations of motion and constraints (in conformal gauge) take the standard form 

- X"^' + r^^(X^X'^ - X'PX"") = 0, (2) 
g^^X^^X"^ = (7^.(X^X- + X'^X"^) = 0, (3) 

where dot ' and prime ' denote derivatives with respect to string worldsheet coordinates r 
and a. Using Schwarzschild coordinates for the background metric 

ds^ = -il )dt^+(l dr^ + r^ dQ^, (4) 



and parameterizing the circular string by ansatz 

t = t{T), r = r(r), 9 = 9{r), 4> = a, (5) 
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one finds that the string equations of motion (Q) and constraints (Q) lead to the following 
system of ordinary differential equations 

*-^(i-^r. (a) 

f = (r - 3M)e^ -{r-M) sin^ 9, (7) 
■■ 2 • 

e = — - sin ^ cos 6, (8) 
r 

supplemented by the constraint 

f2 + (r^ - 2Mr) 0^ + sin^ ^) = E^. (9) 

The integration constant E, appearing here, will play the role of external control parameter 
( "order parameter" ) for the system. It is equal to the total conserved energy of the string 
divided by 27r/i. In the absence of first integrals other than (|), we are thus deahng with a 
three-dimensional phase space. Notice also that the worldsheet time r is not equal to the 
proper time. The worldsheet time r is however finite for a string falling into the black hole. 
Eqs. (1^-0) are equivalent to the following Hamiltonian system 

with the constraint H = 0. This is very similar to the Hamiltonian describing zero angular 
momentum photons in the Schwarzschild background [l^, except for the last term in the 



potential which is due to the string tension. This is precisely the term that leads to non- 
integrability and, as we shall see, chaos. 

The Hamiltonian system (|l^) has been previously considered and solved in the equatorial 



plane = 7r/2) |T3|, and some trajectories in the general case were computed numerically in 
111. In this paper, we make a more complete analysis of the dynamics associated with the 
Hamiltonian (]T0|). In particular, we shall demonstrate the presence of chaos in this simple 
system. 

Due to the non-integrable nature of the dynamical system the analysis of the 

string evolution was done numerically. The string trajectories for various initial conditions 
were obtained by integrating Eqs. (|^) and using the fifth order embedded Runge-Kutta 
method with adaptive step size control [0. The constraint was used to independently 
check for numerical precision. 

It is easy to see that there are three possible asymptotic outcomes of the string dynamics. 
The string can either fly by the black hole and escape to (r, 6) = (oo, vr), or the string can 
backscatter and escape to {r,9) = (oo,0), or the string can be captured by the black hole 
r < 2M. Some examples of string trajectories illustrating these outcomes are shown in 
Fig. 1^, where we plot the string radius R = r sin 6 (vertically) as a function of Z = r cos 9 
(horizontally, string comes in from the left). In the three examples shown, the string is 
initially collapsed at a position a few horizon radii outside of the black hole. It then expands 
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and propagates to the right — towards the black hole. In Fig. ||a, the string passes the black 
hole but then returns and is captured, r < 2M. In Fig. ^d, the string passes the black hole 
and escapes to the right, {r,9) = (oo,7r). In Fig. Qc, the string passes the black hole, then 
returns and finally escapes to the left, (r, 9) = (oo, 0). 

Besides the solutions with these three asymptotic outcomes, there is an infinite set of 
unstable periodic orbits, which separate the solutions with different asymptotic outcomes 
in the phase space of all solutions. Typical examples of periodic orbits are shown in Fig. |^. 
Besides the ones shown, there are also periodic orbits when the string starts a little further 
away from the black hole and oscillates a number of times before reaching it. There are no 
stable periodic orbits in our system. 

To get a better understanding of the string dynamics, we now consider a two-dimensional 
slice of the four-dimensional space of initial conditions in more detail. It is most convenient 
to fix the constant "energy" E and then to impose one more relation between initial values 
of (r, 6, r, 6) at T = 0. Following the procedure of the basin-boundary method for chaotic 
scattering 0, we then color this two-dimensional slice of initial conditions according to 
the three different asymptotic outcomes mentioned above. Due to numerical reasons, the 
determination of the two asymptotic outcomes corresponding to escape is done not at infinity 
but at some large, but finite, value of r. While in general the string can pass this cutoff 
radius in one direction but change its mind later, for sufficiently large r, this procedure will 
only lead to a wrong color for very few points, as follows from the asymptotic behavior of the 
potential in the Hamiltonian (0) (this should be contrasted with the case of the mixmaster 
universe [|16|, where all the trajectories eventually bounce back). 

In Fig. 1^ we show examples of this coloring procedure applied to the phase space of string 
solutions. Fig. ^ and Fig. ^3 show two different two-dimensional slices of the space of initial 
conditions. Fig. ^ and Fig. ^ show magnifications of the regions of Fig. ^ and Fig. ^ 
indicated in the previous figures. The boundaries between different colors in these figures 
correspond to the unstable periodic orbits. The magnifications of these so-called basin- 
boundaries reveal a fractal structure, see Figs. ^b-d. This provides the coordinate-invariant 
indication that the dynamics is in fact chaotic [0. 

To get a quantitative measure we can also determine the fractal dimension. Consider 
for instance Fig. Qd, which corresponds to the two-dimensional slice in the space of initial 
conditions, given by 

E = 14.0M, ^(r cos ^) = at r = 0. (11) 
ar 

The box-counting dimension is defined by 

B = l.m!^#^. (12) 

^-0 ln(l/e) ' ^ ^ 

where N{e) is the number of squares of side length e needed to cover the basin-boundary. A 
square should be counted only if it contains at least two different colors, otherwise it is not 
part of the basin-boundary. The result of this counting is shown as a plot of In A^(e) versus 
ln(l/e) in Fig. ^. The straight line is a least-square fit to the data points and shows that the 
dependence is indeed linear over the very wide range of resolution. The fractal dimension 
of Fig. Hb, calculated from the highest-resolution grid of 4000 x 3200 points, is 
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D = 1.65 ±0.03. 



(13) 



The non-integer value shows in a coordinate-invariant way |T3 that the basin-boundary is 
indeed a fractal; a so-called strange repellor. The error in the result is due to the finite size 
of the numerical grid and slow convergence rate of definition ([T^) , and can be estimated by 
examining D calculated at lower resolutions. 

It is interesting to also examine the fractal dimension as a function of "energy" E. We 
considered the slice (|TT]) in the energy range E E [0, lO^M] and computed the fractal di- 
mension (|12|) of the basin-boundary (corresponding to Fig. ^3 and its analog for different 
energies) as a function of energy E. The result is shown in Fig. where we plot the frac- 
tal dimension D as a function of ln(ii^/M). For low energies, E ^ 4.37M, the dynamics 
is completely regular in the sense that there is only one asymptotic outcome — namely 
capture. For slightly larger energies, with values of E between 4.37M and 5.67M, escape 
becomes possible but the different asymptotic outcomes are connected by "regular transi- 
tions" , corresponding to one-dimensional curves in the space of initial conditions (and thus 
corresponding to fractal dimension 1). However, at energy E ~ 5.67M, the picture changes 
dramatically from regular to highly chaotic dynamics. In a very narrow energy range, the 
fractal dimension changes abruptly from 1 to approximately 1.6, whereafter it increases 
slightly for higher energies. It must be stressed that the numerical values given above for 
energies separating the different "phases" are slice-dependent. The phenomenon explained 
and demonstrated by Fig. ^ is however generic. 

In conclusion, we have shown in a coordinate invariant way that the axially symmetric 
system of a circular test-string in the Schwarzschild black hole background is chaotic. Chaos 
sets in abruptly at a certain "critical" value of the external control parameter, which is 
related to the conserved energy of the string. Below this critical energy, the dynamics is 
regular. Above the critical energy, the dynamics is highly irregular and chaotic. 

In all fairness, it should be mentioned that we have neglected some physical effects, 
including the backreaction and gravitational radiation of the cosmic string. However, the 
system considered here represents the simplest and most symmetric example of string dy- 
namics in black hole spacetimes and therefore suggests quite generally that string dynamics 
in black hole spacetimes is chaotic (notice however that the problem of stationary strings in 
black hole spacetimes is completely integrable |jl8|). It would thus be interesting to consider 
other dynamical string configurations as well, and we notice that there has recently been 



some interest in the scattering and capture of open strings by a black hole; see Refs. p!9|j20 
and references given therein. 
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FIGURES 



String Loop 



Black Hole 



FIG. 1. Schematic illustration of string loop approaching black hole. 





(b) 




FIG. 2. Possible asymptotic outcomes of the string evolution: (a) capture, (b) escape, and (c) 
escape with backscatter. Axes are (r cos ^, r sin 0), and the thick circle represents the event horizon 
of the black hole. 





(b) (c) (d) 

FIG. 3. Unstable periodic orbits. 
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FIG. 4. Basin-boundaries, plotted for E = 14. OM: (a) slice 9 = with r G [2.0M, 29. 2M] on 
horizontal axis, and — r G [0, i?] on vertical axis; (b) slice I;: (r cos 61) = 0, with r G [2.0M,27.5M] 
on horizontal axis, and 6 G [0, 7r/2] on vertical axis; (c) and (d) show fractal detail of figures (b) 
and (c) respectively. 
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FIG. 5. Number of fractal points as a function of resolution, calculated for slice in Fig. ^3. 
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FIG. 6. Fractal dimension as a function of energy calculated for slice in Fig. ^3. Radial size 
of the box is scaled with energy as r € [2.0M, (2.0 + f|;§ ^)M] to keep the choice of the region 
consistent. 
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